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ABSTRACT 
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ground of arbitrary dimension. This propagator is de Sitter invariant and 
possesses the proper flat spacetime and massless limits. Moreover, the re- 
tarded Green's function inferred from it produces the correct classical re- 
sponse to a test source. Our result is expressed in a tensor basis which is 
convenient for performing quantum field theory computations using dimen- 
sional regularization. 
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1 Introduction 



Vector fields described by the Proca Lagrangian [T]: 

C = -\F^F pa gWg v °^-g - l -m 2 A,A v g^ V=£ , (1) 

have great physical relevance because they arise from spontaneous symmetry 
breaking in gauge theories [2] . Although they have been much studied in flat 
space, little attention has been given to other geometries [3]. In this paper we 
shall construct the Feynman propagator of the Proca field described by (pQ) 
on a D-dimensional de Sitter background. Our motivation was to facilitate 
the dimensional regularization computations in de Sitter background of the 
effective potential [I] and of loop scalar bilinear expectation values [5] in the 
theory of scalar quantum electrodynamics; taking the mass to zero gives the 
photon propagator in Lorentz gauge and that gauge turns out to effect great 
simplifications in perturbative calculations. 

One might think the massive vector propagator on de Sitter background 
was found some time ago in the classic work by Allen and Jacobson [6]. 
However, their result diverges in the massless limit, and fails to agree with 
the known flat space limit. The problem turns out to be a term missing from 
the right hand side of their propagator equation which only affects this one 
case of the many they considered. When the missing term is included the 
resulting solution is almost the same as what Allen and Jacobson got, but it 
remains finite in the massless limit and agrees with the known result in the 
flat space limit. 

Section 2 derives the correct propagator equation. In Section 3 we in- 
troduce a way of parameterizing the propagator that greatly facilitates per- 
turbative computations. The propagator equation is solved in Section 4. 
Section 5 gives four useful limits: the flat space limit, the coincidence limit, 
the massless limit, and the limit of D = 4 spacetime dimensions. Section 6 
checks that the massless, retarded Green's function correctly reproduces the 
response to a point charge on the open coordinate sub-manifold. 

2 The Propagator Equation 

The Proca equations associated with ([1]) are: 

d u (V=g-g up g^F pa )-m 2 ^g-g^A u = . (2) 
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For non-zero mass, the Lorentz condition: 

mVO d^y/^gr Aj) = . (3) 

is a direct consequence of the antisymmetry of the field strength tensor: 

Fp, = d^A u -d u A^ . (4) 

We shall continue to impose this condition even for m 2 = 0; in the massless 
case (pTJ) is gauge invariant and fl3J is a gauge condition. Using Q and 
employing covariant notation allows us to express the Proca equations in the 
form: 







V? = y/=g(U 



2; v 

m o„ 



(5) 



where represents the kinetic operator, □ ^ stands for the standard vector 
covariant d'Alembertian operator: 



□ v A — A >p 
(l v — sifi p 



(6) 



and R^ u is the Ricci tensor. 



The Feynman propagator should be the vacuum expectation value of the 
time-ordered product of two vector fields: 



[x: x 



n 



T(A p (x)A a (x') 



n 



(7) 



Assuming that the vacuum state is de Sitter invariant, the propagator must 
be a bi-vector under de Sitter transformations, transforming according to 
x M on the index p and according to x' M on the index o. It would be simple 
enough to construct the propagator by working out the free field expansion 
for the vector potential and then taking the expectation value indicated in 
(J7j) . We will instead follow Allen and Jacobson [B] in searching for a de Sitter 
invariant solution to the differential equations the propagator must obey. If 
this solution is unique, it presumably agrees with the result of the canonical 
operator construction. 

The propagator (JjTJ) must obviously obey the Lorentz condition fl3} on 
both of its arguments: 







-g(x)g^(x)i 



d;{^-g(x>)g^(x')z{pA a ](x;x')} . 



')} 



(8) 
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It must also invert, up to the factor of —i, the kinetic operator T> " on the 
appropriate subspace of transverse functions. This automatically implies 
that: 



V.fi 



P A U 



{x;x') ^ ig^ v 8 D {x - x') , (9) 

because the quantity on the left hand side is transverse on x'^ whereas the 
quantity on the right hand side is not. The correct right hand side must be 
proportional to the transverse projection operator: 

V? i[ p A v ]{x;x') = ig^5 D (x-x') + yj -g{x) d' v iA(x; x') , (10) 

where iA( the propagator of a massless, minimally coupled scalar 

field in de Sitter spacetime. Transversality ([8]) on x'^ follows trivially using 
the equation for iA{x\ x')\ 

y^g(x)g^(x) d u iA(x;x') } = i5 D (x-x') . (11) 

It should be noted that the correct vector propagator equation (jTUI) is 
nothing but the extension to curved spacetimes of the familiar flat spacetime 
relation: 

[x;x') = i(rj, v -^] 5 D {x-x>) , (12) 



(d 2 - 



m 2 



The extra term on the right hand side of this relation, and of ffTUl) . is clearly 
essential for transversality. However, it is just as clearly non-zero away from 
coincidence, that is, for x' M 7^ x M . The technique of Allen and Jacobson was 
to solve the various propagator equations away from coincidence, with zero 
on the right hand side. They then adjusted the coefficients of the linearly 
independent, homogeneous solutions to enforce two conditions [6]: 

• Analyticity at the antipodal point, x ,fl = and 

• The proper singularity structure at coincidence to result in a delta 
function on the right hand side. 

This technique is completely valid when the right hand side contains only 
a delta function, as is the case when a gauge fixing term is added to the 
action. However, it fails for an exact gauge such as ([3]), for which the right 
hand side of the propagator equation must project onto the subspace of 
functions respecting the gauge condition. In that case the right hand side of 
the propagator equation is generally non-zero, even away from coincidence, 
and the propagator must include an inhomogeneous solution to reproduce 
this non-zero term. That is the only thing Allen and Jacobson missed. 
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3 The y-Basis 



Although much of our analysis will be independent of the coordinate sys- 
tem, one should think in terms of conformally flat coordinates for which the 
invariant line element is: 

ds 2 = a 2 (rj) ( — drj 2 + dx ■ dx^j , (13) 
where a(r]) and H are the scale factor and Hubble parameter respectively: 

a(v) = -^ , H 2 ^^— i A>0 . (14) 

Although any scalar function of x^ and x'^ could be expressed in terms 
of the geodesic length l{x ; x'), it is often simpler to use the following length 
function: 

y(x;x') = H 2 a(r]) a(r)') Ax 2 (x ; x') , (15) 
where we have defined: @ 



Ax = y||x — x'|| 2 — (\r] — if\ — is) 2 . (16) 
The relation between y(x ; x') and £(x ; x') is: 

y(x;x') = Asm 2 (^Hiix;^ . (17) 

One can see from (TTTT) that the length function between a point and its 
antipodal is: 

y{x-x) = 4 . (18) 

It has long been realized that any de Sitter invariant bi-vector can be 
expressed in terms of the geodesic length £(x;x'), the gradients of £(x;x') 
with respect to x^ and x ,fl , and the parallel transport matrix \pg u ] x') 
P,[7j. However, this is not convenient for theories such as scalar quantum 
electrodynamics and gravity which possess derivative interactions. The more 
convenient tensor basis consists simply of de Sitter invariant derivatives of 
y(x; x'). Using this basis, derivatives of functions of y(x; x') immediately pro- 
duce basis tensors, and one can dispense with the tedious step of expressing 



1 In (fT4|) A is the cosmological constant. 

2 The infinitesimal imaginary part of (116|) enforces Fcynman boundary conditions. 
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these derivatives in terms of gradients of £(x; x') - using relation (ll?p - and 
the parallel transport matrix - using the identity [8]: 

d 2 y 1 dy dy 



v9v 



[x\ x') 



1 



(19) 



2H 2 dx^ dx' v 2H 2 (A - y) dx» dx' v 

It is simple to work out the contraction of any two basis tensors in con- 
formal coordinates. Because these contractions are de Sitter invariant, the 
relations one obtains in this fashion must then apply in any coordinate sys- 
tem. The contractions we require are [8j: 

dy dy 



g^(x 



dx^ dx v 
dy d 2 y 

dx^ dx v dx' p 
d 2 y d 2 y 

dx p dx' p dx v dx' u 



H\Ay-y' 
H 2 {2-y) 



dy 

dx'p 



AH 4 g pa (x') - H' 



dy dy 
dx' p dx" 7 



We will also need covariant derivatives of the basis tensors [SJ: 



D dy 



d 2 y 



Dxp dx v 
v d 2 y 



dx^dx h 



dy 

dxP 



H 2 {2-y) 9ll] 



x 



-H 2 



d 2 y 



(20) 
(21) 
(22) 

(23) 
(24) 



M dx u dx'p " dxvdx'p 

Each of these identities ( 12011241) applies as well when x^ and x'p are inter- 
changed. 

One can express any de Sitter invariant vector propagator as follows in 
the y-basis [SJ: 



[ - 0C J 



B{y) 



d 2 y 



+ C(y) 



dy dy 



(25) 



dxP dx' v dxP dx' v 

Substituting this representation into the transversality relation ([8]) and mak- 
ing use of identities (I20H24I) gives: 

-DB(y) + {2-y)B'{y) + {D+l){2-y)C(y) + (iy-y 2 )C'{y) = . (26) 

This equation implies that the functions B(y) and C{y) can be expressed in 
terms of a single function j(y): 



B{y) 
C(y) 



l 



4(D - 1)H 2 L * 'dy 

1 r d 

\(?-y)»r-( D 



{Ay-y 2 )^- - {D-l)(2-y)h(y) , (27) 



A(D - 1)H 2 



dy 
5 



i{y) ■ 



(28) 



Relations (I25H28I) are independent of the propagator equation (TTUj) . In fact 
these formulae are simply the translation into the ?/-basis of results previously 
obtained by Allen and Jacobson [6]. 

We close this section by expressing the propagator equation f]T0]) in the 
y-basis, beginning with the right hand side. It has long been recognized that 
equation ffTT]) for the scalar propagator has no de Sitter invariant solution 
[9]. However, the mixed derivative <9 M c^iA(x; x') which appears in the vector 
propagator equation (flQj) is nonetheless de Sitter invariant, up to a delta 
function term. H If we elect to maintain homogeneity and isotropy, the 
scalar propagator takes the form jlOj : 



iA(x; x') = A(y) + kin a(r]) a(rj' 



(29) 



Therefore, the extra term on the right hand side of ( fTOj) is: 

— x') 
dy dy 



-g(x) d^dl iA(x;x') 



5y v i5 D {x-x') 



A"(y) 



dx^ dx' v 



+ A'(y) 



d 2 y 



dxv- dx' v 



• (30) 



The de Sitter breaking logarithms of ([29]) drop out in the mixed derivative, 
and we need not worry about the delta function terms if we follow Allen and 
Jacobson in solving the equation away from coincidence. 

Before going on to the vector propagator it is worth pointing out a few 
more things about the scalar propagator. First, its equation (TIT]) implies the 
following relation for the function A(y) for y ^ [10J: 

y 2 )A"{y) + D(2 



(Ay-y 2 )A"(y) + D(2-y)A'(y) = (D - 1) k 
Second, the de Sitter invariant function A(y) is [11] : 



(31) 



A(y) 



H 



D~2 



(47T 



2 



+ E 

ra=0 



:f-i) 



r(f + l + n) 



+ 



r(f + i) 
(f-2) 

7TCOt(^ J D) 




r(f; 



(«)- 



'r+2 



r(n + 2) f - 2 - n 



+ 



T(D-l+n) (I)" 



r(f + n) 



n 



,(32) 



3 For example, it is the coincidence limit of this mixed derivative that occurs in the de 
Sitter invariant expectation value of the scalar stress tensor [9]. 
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and the constant k has the value [10J: 

H D-2 pp _ -Q 



fc = 



(4tt 



r(f: 



(33) 



Finally, the derivative of A(y) can be expressed in terms of a hypergeometric 
function [T2] : 



1 r(D)T(l) / D y 

A{y) = TcfTiy 2Fi r >1 'T +1 ' 1 "4 



(34) 



We are now ready to tackle the left hand side of the propagator equation 
fTTOj) . By substituting (125]) and (130]) . and again making use of (I2QH24|) . we 
obtain two differential equations involving B(y) and C(y) away from the 
singular point y = 0. The term proportional to (d 2 y /dx dx') implies: 



m 



1 



(4y-y 2 )B" + D(2-y)B'- \ D + — J B + 2(2- y)C = —^A' , (35) 
while the term proportional to (dy/dx) (dy/dx') leads to: 
(Ay-y 2 )C" + (D + A)(2-y)C- ( 2 D + ^\c-2B' = ± 



A" . (36) 



Now B(y) and C(y) are not independent variables; they are related to one 
another by the transversality constraint ( 126]) . Using ( |26]) and ( |3T]) . it is 
straightforward to see that the following combination of equations (135]) and 
(136]) gives a tautology of the form = 0: 

d d 
(2 - y) — (eqrm> + (4y - y 2 ) T-(egri36]) - D (eqr^ 
ay ay 

+ (D + l){2 - y) {eqrMI ■ (37) 

The independent combination of (135]) and (136]) can be obtained by invert- 
ing relations (1271) and (128]) to express j(y) as: 



7(2/) = -H 2 
Its derivatives equal: 



{2-y)B(y) + {Ay-y 2 )C{y) 



7 



7 



(2 - y)# -B + (4y- y 2 )C + 2(2 - y)C 



(38) 



(39) 



(2 - y)B" - 25' + (4y - y 2 )C" + 4(2 - y)C - 2C . (40) 
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Now we make use of (13"5]1 . (1551) and (13"T1) to infer: 

{Ay-y 2 )i\y) + (£> + 2) (2 - y) 7 '(y) 

= -(4y-y 2 )A"(y) - (2-y)A'(y) 



m 



i{y) 



(D-l) \(2-y)A'{y) - k 



(41) 
(42) 



This is the fundamental expression of the propagator equation (TlOl) in terms 
of the y-basis. 



4 The Solution 

The general solution of equation ( 1421) takes the form: 

l{y) = Ci7i(y) + c 2 72(y) + J P (y) , (43) 

where 71 (y) and 72(1/) are the linearly independent solutions of the homo- 
geneous version of equation (T42l) and j p (y) a particular solution of the full, 
inhomogeneous equation. From (T3TT) it is straightforward to see that a par- 
ticular solution is: 

lp {y) = — 2 (D - 1) [-(2 - y)A\y) - k] . (44) 

Using relation ( 1341) and a few simple facts about hypergeometric functions 
[12], we can write 7p (y) explicitly as: 

H 2 H D - 2 T(D) _ / D y 



7p(y) = -^77T1 ^ ^ 2*1 U>-l,2,- + l;l-f . (45) 



(4tt)t 2T(f + 1) * 1 V ' ' 2 4 



The two homogeneous solutions can be similarly expressed as hypergeometric 
functions IT2l: 



^to) - ^(^ + <^-, § + i;i-f) • («) 

= ,*i (E±± +v ,£±±- v ,E + 1 ;l) , (47, 



where we have defined: 



D-3\ 2 m? 
2 ) ~~ IP 



(48) 
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We must now determine the two arbitrary constants C\ and c 2 present in 
by enforcing the proper behaviour at y = and at y = 4. 

The desired behavior at y — 4 - that is, at the antipodal point - is for 
the solution to be analytic. This is most easily examined by employing the 
appropriate expansion of the hypergeometric function 



2 Fx (a , , 6 ; z) 



E 

n=0 



T(a + n)T(0 + n)T(5) z^_ 
T(a)T(0)T(5 + n) n[ 



a0 a(a + 1)0(0+1) ^ 3 
+ 5 + 5(5 + 1) 2! + U[ J 



(49) 
(50) 



It is apparent that r ) p (y) - given by (|4"5]) - and 71 (y) - given by (146]) - are 
analytic at y = 4 while 72(1/) - given by ( l4"Tj) - is not. We must, therefore, 
set: 

c 2 = =^ 7(2/) = c X 7a(2/) + 7 P (y) . (51) 
The desired behavior at y — - that is, at coincidence - is for the solution 



to be ~ y 2 . This is achieved by enforcing the following condition on (l5Tj) : 



7 



Cl 



7i 



(*)T 



2 + 



7 P 



\ J/ / 







(52) 



where the subscript on the brackets in f)52p indicates the term in the Laurent 
expansion of the bracketed quantity which is proportional to the subscript. 
To determine this part we use one of the transformation formulae of hyper- 
geometric functions 



2 F 1 (a , , 5 ; z) 



T(5)T(5-a- 0) 



| T(5)T(a + 0-5) (i 



T(5-a)T(5-0) 

\5-a-P 



2F1 (a , , a + - 5 + 1 ; 1 



T(a)T(0) 

x 2 Fi(5-a,5-0,5-a-0 + l; 1 



(53) 



Applying ( 1531) to our particular values: 



D + l 



a 



+ v 



(3 



D + l 



— v 



2 11 2 

gives for the relevant part of our hypergeometric function 



/D + l D + l 

l F l o H v > o 



v , h 1 ; 1 

'2 ' 4 



(54) 
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r(f + i)rffi 



D_ 

2 



(55) 



Moreover, using (jMj) we also conclude 
A' 



1 £^ £ ^ 
4 ( 47r )f V vj/, 



(56) 



Inserting (I55H56P into (152]) reveals that the latter can only be satisfied if: 
D-l H 2 H D - 2 r(^±l + „) T(2±l - v ) 



Ci 



(47T) 



r(f + i) 



(57) 



Having determined both constants C\ and C2, we arrive at the unique 
solution: 

1 H 2 H D ~ 2 



D 



+ 



(Any 

D+l 



2F Jo-l,2° 



D 
D+l 



r(f + i) 



^-1 



r(f + i) 



(58) 



x 2J Pi 



2 

£ + 1 



7/ . 



£ + 1 D y 
— : V , — + 1:1 — - 



where the parameter v is defined in (I48p . The second hypergeometric function 
in (|58|) represents the solution obtained by Allen and Jacobson [BJ; the first 
hypergeometric function is our contribution. It is often useful to have the 
Laurent expansion of (|58|) : 



i{y) 



D-l H 2 H°- 2 



+ E 

n=0 



(47T) 



m 2 n 
# 2 1 2 



1) 



(n + 1) 



T(n + D- 1) 
T(n + f + 1) 



+ 



r( 



+ (n 

i)r(2-f) 



~2 



r(| + 1/) r(| - 1/) 



E 



n=0 
5 
2 



+ 3) 
T(n 



r(n + f + 1) 



rc-f+2 



I> + 3) 

iai + z/ )r(n + ^±i 



r(n+f + l)T(n + l) 



r(n+| + z/)r(n + |-^) 
I> + 3)I> - f + 3) 



(59) 
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5 Various Limits 



In this section we derive interesting limits for the function 7(2/) whose al- 
ternate expressions, (!58|) and (|59|) . have just been given. Where appropriate 
we also derive limits for the full propagator, whose dependence upon 7(7/) is 
restated here for convenience: 



[X\ X 



B{y) 



d 2 y 



+ C(y) 



dy dy 



B{y) 
C(y) 



4(D-1)H 2 
1 

A(D- 1)H 2 



dx^ dx' u dx^ dx 1 v 

2\ 



i±y-y 2 )i{y) - {D - 1)(2 - y)j(y) 



(2-y) 1 \y) - (D-l)j(y) 



(60) 
(61) 

(62) 



5.1 Flat Spacetime Limit 



One takes the flat space limit by letting H approach zero with the physical 
time t fixed, rather than the conformal time 77: 



Flat 



H 



& t fixed 
1 



-Ht 



H 



+ t 



(63) 
(64) 



We will express the flat space coordinate using the symbol X M , meaning 



X° = t and X 1 



x 



From the definition (|15l) of y(x; x'), and relation f[6"4"|) . 



we see that the flat space limiting forms for y(x; x') and its derivatives are: 
Flat =■ //(./•:./•') — > H 2 AX 2 , (65) 

-> -AH^AX^AX, , (66) 



dy dy 



dx 11 dx' u 
d 2 y 



(67) 



dx^ dx' v 

Here rj^ is the Minkowski metric and the Minkowski coordinate interval is: 



AX 



V,u(X-X' 



\x-x 



'112 



(\t-t'\ -ie)' 



(68) 
(69) 



We also need asymptotic forms for functions of the parameter v in the 
limit of small H with the mass m fixed. From expression (HHj) we see: 



im 



(70) 
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Now use the Stirling approximation to obtain the asymptotic form of the 
ratios of Gamma functions which occur in the Laurent expansion (1591) : 



T{K+\±v 



1 I im n 

2 H i 



p/1 ± vm\ 



±- 



im\K 



H 



(71) 



By using (1711) alternatively for = n + y and X = n + 2, we compute the 
flat space limit of (1591) : 



Km 7(2/) 

flat 



D - 1 ,£> x / £>x m D ~ 2 
2 v 2" J v 2 ~~ ~2> (47r)f 



00 

xE 

ra=0 



;> 2 Ai 2 ) n ^ +1 



;±m 2 AX 2 )™ 



r(n + 2- f)(n + l)! r(n + 1 + f )n! 



(72) 



We shall employ a tilde to denote the flat space limit of any quantity; for 
instance: 

7(AX 2 ) = IWy) . (73) 

/<at 

It is straightforward to check that expression (1721) obeys the equation: 



4AX 2 f(AX 2 ) + 2(^ + 2)7 (AX 2 ) - m 2 7(AP 



D-i r 



The analogous quantity 7 for Allen and Jacobson is: 



7a.,(AX 2 ) = 7(AX 5 



D — 1 m 



D-2 



7TT AX^ 



(74) 



t (±m 2 AX 2 )f 



2 (4vr 

This function obeys the homogeneous analogue of (1741) : 
4 AX 2 7L(AX 2 ) + 2P + 2) 7L(AX 2 ) - m 2 7 AJ (AX 2 







(75) 



(76) 



The flat space limit of the full propagator additionally requires expressions 
(E0H62D and the limiting forms {66117}: 



(x: x 



Vfl „ 7(AX 2 ) + 



D - 1 



T]^ AX 2 — AX M AX U J 7 ; ( AX 2 ) . (77) 
12 



Using relation (J73|) it is easy to check that (I77|) obeys the correct equation: 



(«9 2 - m 2 



x; a; 



d 2 



5 D (x — x') . 



(78) 



The right hand side of (1731) results in the non-local term on the right hand 
side of (17H|) which is non-zero away from coincidence. The flat space limit 
of the Allen- Jacobson propagator is obtained by replacing 7 by 7 AJ in (177|) . 
In contrast to the correct propagator, one can use relation (J75|) to show that 
the flat space limit of the Allen- Jacobson propagator obeys: 



{d 2 



m 2 ) i 



A AJ 



[x — X 



(79) 



The extra term in ( 1751) has promoted the singularity from a delta function 
to second derivatives of a delta function. 



5.2 Coincidence Limit 

Obviously "coincidence" means = This implies: 

• All positive powers of y vanish; 

• All D-dependent powers of y vanish when using dimensional regular- 
ization; and 

• The tensors of the y-basis have the following limits: 

dy dy d 2 y 



lim 

x'-*x ox^ OX 1 ' 







lim 

x'-+x ox^ox' 1 



2H 2 9lMU . (80) 



The first and second rules imply that only the n = terms from the two 
infinite sums of (1591) survive: 

D-IH 2 H D - 2 



7(0) 



2 m 2 (4vr 

r(f + 1) 



(81) 



r( 



' "2' 



- 1/ 



The first term inside the square brackets derives from our correction. The 
coincidence limit of the full propagator requires the third rule, in addition to 
relations (160|). (IBTD and (1621): 



lim i 

X 1 — >x 



7(0) g^u 



(82) 
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5.3 Massless Limit 



This is the most important limit because it gives the Lorentz gauge propaga- 
tor. We therefore distinguish it with the convention that an overline indicates 
the massless limit: 



lim n 7(2/) = l(y) 

m— >0 



lim i 

m— >0 



[x\ x') 



(83) 



To take the massless limit we first expand the parameter v (148]) about its 
value for m = 0: 



D-3 



v 



1 m _ ( m 

T7K + O — 



D-3 H 2 



i 

H 1 



D-3 



+ Av . 



(84) 



We then substitute into form (|59|) for 7(2/) and expand the ^-dependent 
Gamma functions using the polygamma function [12] : 



m = ^ln[r(z)] = ip(z — 1) + ~~y 



(85) 



The result is: 



i(y) 



D-l H 



D-2 



2{D - 3) (4 



7T 2 



+ E 

n=0 



(n + l)V(n + D- 1) 
r(n + f + 1) 



^(2 - y) - V(§ - 1) + + D - 1) - ^(n + 2) 



E 

n=0 



n 



D 



+ 3)T(n + f + 1) 
T(n + 3) 



(86) 



L> D D D 

^(2 - -) - V(- - 1) + ^(n + - + 1) - *P(n - - + 4) 



n-f+2 



Note that 7(7/) would diverge without the extra terms we have added to the 
corrected solution ( 1581) . The result for the massless limit of the propagator 
comes from simply replacing 7(2/) everywhere by j(y) in expressions (|60|) . 
(ED and (E21). 
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Additionally taking the flat space limit gives: 

D-irfS-1) 1 



Km 7(0) 

flat 



■ D N 



limi 

flat 



[x: x 



Air 2 



47r f 



+ 



D-2 ' AX 2 
Additionally taking the coincidence limit gives: 
D - 1 H°- 2 T(D - 1) 



7(0) 



2(D-3) (4vr)# r(f + 1) 
4>(D-l) - ^(2) 



1 



+ ^(2-§) 



(87) 
(88) 



The coincidence limit of the massless limit of the propagator comes 
simply replacing 7(0) by 7(0) in expression 



(89) 
from 



5.4 D = 4 Limit 



When using dimensional regularization in perturbative computations involv- 
ing scalar QED, we can take D = A on any part of the result which makes a 
finite contribution. This means one typically only needs to keep D arbitrary 
in a few of the terms of expression (1861) for 7(2/). We have found the following 
expansions useful pfl [5] : 

7(2/) = 



and: 



7(l/) 




D 



-r(- 

V2 



L>-2 



(Air) 



D-2 



HI) 



y\2 
a) 



+ 



O(D-A) 



(91) 
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6 Response to a Point Source 

We begin by extracting a pure gauge term from our propagator (!60|) using 
relations (EHES and (SI: 



;x; a; 



1 d j dy 



+ 



(Ay - y 2 ) i(y) 
+ (D-l)(2-y)j(y) 

4{D-l)H>dxH*» I* 

+ (D + 2)(2-y) 1 '(y) - 2{D-l) 1 {y) 
1 d j dy 



dy dy 



(92) 



A(D - 1)# 2 <9x^ [&r /! 



1 



4# 2 



1 



m 



(4y - y 2 ) j'(y) 

+ (D-l)(2-y) 1 (y) 
dy dy 



7 (l/) + (2 - y)A'(y) - k 



dx^ dx' 



(93) 



D — 1 H 2 

Now recall that the retarded Green's function derives from the imaginary 
part of the propagator: 

/?f](x;x') = -29(ri-ri , )lm{i[pA v ](x;x')} . (94) 

It only makes sense to extract the pure gauge term in the massless limit, for 
which the theory is gauge invariant. When we also take D = 4, we conclude: 

lim (2-y)A'(y) = -|1 {-A + 2) . (95) 



* ' 16tt 2 
Hence the retarded Green's function is: 



■ret' 



(x; x') = (pure gauge) — 



29(Ar]) H 2 



AH 2 16tt 2 



Im< — 



dy dy 



(y + ie) 2 j dx^ dx' 1 



(pure gauge) + — (y) 



(96) 

4tt " N< " dxv dx' v ' 
The retarded Green's function determines the field response A^(x) due 
to the presence of a source J v (x): 



j4 / 

a x 



ret 



(x;x') J u {x') 



(97) 
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We shall be concerned with a point source so that: 

r(x) = q fdrx u (r)S\x- X (r)) 



(98) 



where x M ( r ) is the geodesic of the source. Up to a pure gauge term, the point 
source response successively equals: 



d A x' 6{At]) 5'{y) 



dr 



q_ 

An 

— I dr — 

An 

An 

q_ 

An 



dy 
dxt 1 



dr^r v nr)5\x'- X {r)) , 



y(x; X (r)) e( V - X °(r)) 5'(y) 



dy 

dx^ 
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d dy 



dy 
dx^ 



dy dx>- 



J y=0 



For a particle at rest on the conformal origin the length function is: 

x ^t)=t5% =► y(*;x(r)) = 



x 



(v 



(99) 



(100) 



where r represents the conformal time of the source and r\ the conformal time 
of the observation event. The relevant derivatives for our computation are: 



dy_ 

dr/ 



x 1 ^ r 



dy 

dx i 



2x l 
tjt 



Tj T T T) 

leading to the final expressions for the response: 



dy 
dr 



2x 



y=0 



77(77 — x) 



(101) 



A (x) 
Mx) 



An 



x 2 + if — r/x 
r/x (tj — x) 



X 



Air x (77 — x) 
and invariant field strength tensor: 



q x 
An x 3 



(102) 
(103) 

(104) 



This is the correct result in our coordinate system [13] so we conclude that 
our propagator provides the proper response. 
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7 Epilogue 



We have worked out the Feynman propagator for a massive vector field on 
/^-dimensional de Sitter space. This involved repairing a minor error in the 
classic analysis of Allen and Jacobson [6J. They failed to include a term 
which is non-zero even away from coincidence and is required to make the 
right hand side of the propagator equation ([TO]) respect the transversality 
of the left hand side. Without this term, the propagator diverges in the 
massless limit and also fails to agree with the known result in the flat space 
limit. The correction is nonetheless small, and much of the solution is simply 
the transcription to our notation of the work of Allen and Jacobson. 

Notation deserves comment as well. It is traditional to represent de Sit- 
ter invariant bi-vectors such as the propagator with scalar functions of the 
geodesic length, i(x\ x'), times a tensor basis formed from gradients of £(x; x') 
on x^ and and the parallel transport matrix x'). We have instead 

chosen to use a function y(x; x') of the geodesic length which has a simple 
expression (j!5p in the conformal coordinate system most often used for ex- 
plicit computations. We have also taken our tensor basis to consist of the de 
Sitter invariant first and mixed second derivatives of y(x;x'). This greatly 
facilitates computations in theories with derivative interactions such as scalar 
quantum electrodynamics [H [5] and gravity. A few simple identities suffice 
for all contractions (I20H22P and covariant derivatives (12311241) of these tensors. 

Finally, we demonstrated that the retarded Green's function inferred from 
our propagator (in the massless limit) correctly reproduces the known clas- 
sical response to a point charge on the open sub-manifold. This is an issue 
because the peculiar causal structure and linearization instability of de Sitter 
space invalidate some of the familiar features of flat space electrodynamics 
[HI [EE ETUI EH] • There also seems to be a poorly understood obstacle to im- 
posing certain gauges in gravity and electrodynamics [18], [13]. In particular, 
it has been shown that de Sitter- Feynman gauge leads to on-shell singularities 
in the one loop scalar self-mass-squared of scalar quantum electrodynamics 
[S] . It is reassuring that our propagator leads to no such problem [5] [JJJ] . 
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